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We study the possibility to replace the anti-neutrino run of a long baseline neutrino oscillation 
experiment, with anti-neutrinos from muon decay at rest. The low energy of these neutrinos allows 
the use of inverse beta decay for detection in a Gadolinium-doped water Cerenkov detector. We 
show that this approach yields a factor of five times larger anti-neutrino event sample. The resulting 
discovery reaches in 613 , the mass hierarchy and leptonic CP violation are compared with those from 
a conventional superbeam experiment with combined neutrino and anti-neutrino running. We find 
that this approach yields a greatly improved reach for CP violation and 6'i3 while leaving the ability 
to measure the mass hierarchy intact. 

PACS numbers: 14.60.Lm, 14.60.Pq, 29.20.dg, 29.40.Ka 



The fact that our Universe is made out of matter and 
not anti-matter constitutes one of the great open ques- 
tions in physics. The conditions required to obtain a 
baryon asymmetry from a CP symmetric initial state 
were identified by Sakharov more than 50 years ago [1|. 
They are deviation from thermal equilibrium, baryon 
number violation and CP violation. Of these three in- 
gredients, CP violation has proven to be the most elu- 
sive. While some small amount of CP violation exists 
in quark mixing [2|, |3| , it is not sufficient to explain the 
baryon asymmetry of the Universe [J]. Therefore, the 
quest for new sources of CP violation continues to be a 
guiding theme for much of particle physics. Leptonic CP 
violation has been proposed as a mechanism to explain 
the baryon asymmetry of the Universe, so called leptoge- 
nesis [5| , and with the discovery of neutrino oscillation [te] 
leptogenesis has become a very plausible scenario. The 
discovery of CP violation in the lepton sector would con- 
stitute a smoking gun for this mechanism. In addition the 
mere existence of neutrino oscillations implies a neutrino 
mixing matrix in the same fashion as the CKM matrix 
in the quark sector, and with three generations there will 
be at least one CP phase. 

Studying CP violation in neutrino oscillation, requires 
the use of appearance channels, when the neutrino 
changes flavor from production to detection. Since i^r 
are very hard to detect and to produce, naturally, all 
proposed methods focus on using transitions between 
Vf, -f-^- Vf^i and Pg -f-^- i^^j. There are various proposals for 
novel techniques to produce neutrino and anti-neutrino 
beams, either from the decays of muons (neutrino facto- 
ries), or from the /3-decay of short-lived, artificially pro- 
duced radioactive nuclei (/3-beams), see e.g. Q- While 
these concepts may ultimately provide the best sensitiv- 
ity to CP violation, in the near term conventional neu- 
trino and anti- neutrino beams based on the decay of horn 
focused pious are the option of choice (superbeams). All 
beam based approaches to neutrino oscillation are per- 
formed with neutrinos in the energy range from a few 



hundred to many thousand MeV . In this energy range 
charged current reactions provide the dominant detec- 
tion mechanism. At the atmospheric mass squared dif- 
ference these energies imply a source detector distance 
from 100 km up to several 1,000 km, hence these experi- 
ments are known as long baseline experiments. 

In the context of superbeam experiments, a CP viola- 
tion measurement requires data from both j/^ -H> i^e and 
Df^ ^ De- However, the anti-neutrino run poses a number 
of specific challenges: the production rate for tt^, the par- 
ent of v^, is lower than for tt"*", the anti-neutrino charged 
current cross section is lower, the background levels are 
higher^, and the systematic errors are expected to be 
larger. Overall, the event rate for anti-neutrinos is sup- 
pressed by a factor of 2-5, depending on the anti- neutrino 
energy, which is illustrated by table H] 

On the other hand, for energies below about 100 MeV 
the situation is quite the opposite, strongly favoring the 
use of anti-neutrinos. In this energy range the dominant 
process is inverse /3-decay, (IBD), i/g -f p — >■ e"*" -I- n. This 
process provides a very useful delayed coincidence tag 
between the prompt positron and the delayed neutron 
capture, which is the reason why this reaction was used 
by Reines and Cowan to discover the neutrino [Sj. At 
these energies copious sources of anti-neutrinos exist. In 
this letter wc focus on j/^ from /i+ decay at rest from a 
stopped pion beam. 

In a stopped pion source a proton beam of a few GeV 
energy interacts in a target producing both 7r+ and tt~ . 
The pions then are brought to rest in a high-Z beam stop. 
The 7r~, and also the daughter /i^, will capture on the 
high-Z nuclei. The 7r+ will produce the following cascade 



^ This is due to the larger contamination from wrong sign pions. 
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FIG. 1: Discovery reaches at Scr confidence level from left to right for Qy6, 
DAR+HFN for normal (NH) and inverted (IH) true mass hierarchy. 



mass hierarchy and CP violation for 6 years of 



of decays 



resulting in v^ , P^ and z/g , but no 9^ . The D^ can oscillate 
into a z/g and this De can be detected and uniquely iden- 
tified by IBD. This concept has been used for instance in 
the LSND experiment [9|. We will denote this combina- 
tion of anti-neutrino source and IBD detection as muon 
decay at rest (dar). 

The central idea of this letter is to combine a horn fo- 
cused high energy i^^ beam, henceforth denoted as horn 
focused neutrino beam (hfn) with i/^ from a DAR setup 
to study ^13, the mass hierarchy and leptonic CP vio- 
lation. This is in contrast to reference [lOJ, where an 
entirely dar based experiment was considered to mea- 
sure leptonic CP violation. In a pure DAR setup there 
is no information on the mass hierarchy and therefore a 
CP measurement is only feasible assuming that the mass 
hierarchy has been determined by other means. We show 
that the combination of HFN and DAR allows one to ad- 
dress both the mass hierarchy and CP violation at the 
same time. Throughout this letter, we will denote this 
new technique as DAR-f hfn. 

To illustrate the strength of dar+hfn, we will study 
a specific setup, which closely resembles the Fermilab 
DUSEL concept for a long baseline experiment, currently 
known as LBNE. Obviously, similar considerations hold 
for any superbeam experiment. This setup has a total 
running time of 6 years and a 100 kt water Cerenkov de- 
tector. The entire HFN part is very similar to the setup 



described in detail in [11[ , specifically we take the source 



detector distance to be 1300 km and use the same detec- 



tor performance. The beam delivers 6.2 x 10 protons on 
target per year, which for 120 GeV protons corresponds 
roughly to 700 kW of beam power. 

The feasibility of a high intensity DAR source hinges on 
the availability of cost effective, megawatt proton accel- 
erators in the GeV range. We use the same proton source 
parameters as in reference [10| which are 4 x 10^^ of i^e, 
v^ and v^ per flavor per year per accelerator. We use 4 
of these accelerators with a source detector distance of 
20 km. We determined this distance to be the optimal 
baseline for our purposes, by performing a scan of base- 
lines in 1 km increments from 4 — 24 km. We also tested 
that putting all sources at the same distance provides 
the best performance. In order to allow for an efficient 
detection of the IBD signal, the water Cerenkov detec- 
tor will need to be doped with Gadolinium (Gd) Il2| . 
Gadolinium- doping decreases the neutron capture time 
and improves the neutron capture signature. Recent ex- 
perimental tests of this concept [13| indicate a detection 
efficiency of 67%, which we will use for our calculation. 
Due to the strong tt" absorption in the target, the Pe con- 
tamination is very small ~ 4 x 10~^, but fully taken into 
account [9|. The non-beam backgrounds, which mostly 
stem from atmospheric neutrino interactions are taken 
from reference [lO| and scaled to a 6 year run with a 
100 kt detector. The energy resolution for IBD events is 
identical to the one for positrons and is parametrized as 
a{E) = 5%y/E/GeY M- 

The very different duty factors, d, of the HFN source, 
d < 10^"*, and DAR source, d ~ 0.1, allow for concurrent 
operation. As shown in table U the event numbers from 
both sources are comparable and thus the probability of 
finding events from both source happening at the same 



time is approximately given by the ratio of their duty 
factors, which is 0.001. This is equivalent to less than 0.5 
events for the duration of the experiment. Obviously, for 
a superbeam it is not possible to run neutrinos and anti- 
neutrinos simultaneously, since the horn^ either focuses 
7r+ or TT^. 

Note, that the disparate baselines reflect the different 
neutrino and anti-neutrino energies and the ratio of base- 
line to energy (L/E) is nearly the same for hfn and DAR. 
For comparison we also present results for a pure super- 
beam experiment of the same total duration, target mass 
and number of protons on target per year, which consists 
of 3 years of hfn operation followed by 3 years of horn 
focused anti-neutrino (hfa) operation. The performance 
for both periods is taken from reference ll|. This setup 
is denoted by hfa-|-hfn. 

In table U we compare the total event rates for the 6 
year DAR-|-hfn and hfa-|-hfn runs. For the oscillation 
parameters chosen, DAR-|-hfn has twice the statistics in 
the neutrino mode and five times as much statistics with 
a five times better signal to background ratio in the anti- 
neutrino mode. 
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TABLE I: Comparison of the signal and background event 
rates of 6 years running of dar-|-hfn and hfa+hfn. Note, 
that for DAR-f hfn this is 6 years of simultaneous running 
of neutrino and anti-neutrinos, whereas for hfa+hfn tliis is 
3 years each, run consecutively. Oscillation parameters are 
sin^ 2^13 = 0.1 and (5cp = and normal hierarchy. 



For our statistical analysis we use the techniques out- 
lined in [l5| . We bin our DAR data into bins of 1 MeV 
and the HFN/hfa data into bins of 125 MeV. We in- 
clude a 5% systematic on the total number of signal 
events and a (uncorrelated) 5% systematic on the total 
number of background events. All physics sensitivities 
have been computed using GLoBES |16| and are shown 
at 3 (T confidence level. The oscillation parameters are 
sin2 6li2 = 0.3, sm^e23 = 0.5, Am§i = 2.4 x 10-3eV^ 
Am2i — 7.9 X 10~^eV^ and we include a 10% error on 
the atmospheric parameters, a 5% error on the solar pa- 
rameters and a 2% matter density uncertainty. In the fit 
the mass hierarchy is free. 

In figure [1] we show the resulting discovery reaches for 
013, the mass hierarchy and CP violation for dar-I-hfn as 
a function of the true values of sin^ 20i3 and (5cp for both 



the normal and inverted true hierarchy as labeled in the 
legend. Due to the asymmetric nature of this experiment 
and the fact that matter effects are relevant only for the 
HFN run, there is a pronounced effect from changing the 
input (true) mass hierarchy. However, in both cases the 
synergy between the HFN and DAR is apparent. 

In figure [2l we compare the results from dar+hfn 
with hfa-I-hfn. The reaches are given as a fraction of 
i5cp and as a function of the true value of sin 20i3 . In 
panel (a), we show the results for the discovery of the 
6*13 and find that DAR-|-HFN outperforms the superbeam 
experiment HFA-|-HFN for all CP phases and both hier- 
archies by roughly a factor two. The discovery reach for 
the mass hierarchy is shown in panel (b) and here, we see 
that for some values of the CP phase, in particular for 
inverted mass hierarchy, the reach is somewhat smaller 
for DAR+HFN. If at the end of the dar+hfn run, the 
mass hierarchy has not been discovered adding a hfa run 
may be required. Finally, in panel (c) the discovery reach 
for CP violation is shown. DAR+HFN improves the reach 
for small 6*13 by a factor between 3 and 10 depending on 
the mass hierarchy and by about 75% for large 6*13. 

Summarizing, we have demonstrated that a combina- 
tion of low energy i/^ from muon decay at rest with 
high energy Up, from a superbeam aimed at the same 
Gadolinium-doped water Cerenkov detector yields a 
moderately improved reach for ^13 and a vastly improved 
discovery reach for CP violation while only marginally 
affecting the mass hierarchy sensitivity. These improve- 
ments are a direct result of combining an optimized neu- 
trino with an optimized anti-neutrino run. The prac- 
ticality of this proposal, however, depends critically on 
the feasibility of low-cost, high-intensity stopped pion 
sources. 
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^ With a solenoid it is possible to foeus both types of pions, but 
that leaves the problem of Icpton charge identification in the 
detector. 
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FIG. 2: CP fractions for which a discovery at 3 a confidence level is possible as function of sin^20i3. From left to right for 
^13, mass hierarchy and CP violation. The different lines are for normal (NH) and inverted (IH) true mass hierarchies and for 
DAR+HFN and HFA+HFN, respectively. 
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